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We calculate the flavor-singlet contribution to the B —f rj^'^ transition form factors from the 
gluonic content of the Ty''' meson in the large-recoil region using the perturbative QCD approach. 
The formulation for the 7)-r;' mixing in the quark-flavor and singlet-octet schemes is compared, and 
. . . employed to determine the chiral enhancement scales associated with the two-parton twist-3 r;''' 

' meson distribution amplitudes. It is found that the gluonic contribution is negligible in the B —> rj 

I form factors, and reaches few percents in the B ^ r]' ones. Its impact on the accommodation of the 

. measured B — > rf'^K branching ratios in the perturbative QCD and QCD-improved factorization 

(N ■ approaches is elaborated. 
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It is still uncertain whether the flavor-singlet contributions to B meson decays into rj^'^ mesons play an essential 
, role. The flavor-singlet contributions to the B — > ij^'^K branching ratios, including those from the b — > sgg 
' transition [l| , from the spectator scattering 0, Q , and from the weak annihilation, have been analyzed in the 
I QCD-improved factorization (QCDF) approach f4!|. However, at least the piece from the weak annihilation 
■ can not be estimated unambiguously due to the presence of the end-point singularities. The flavor-singlet 
contribution to the B rj^'^ transition form factors from the gluonic content of the ry^'-* meson is also involved 
in the annihilation amplitudes. Though this contribution seems to be most crucial among all the flavor-singlet 
pieces in the B rj'^'^K decays, it has been parameterized and varied arbitrarily The form factors associated 
Q-!. with the decays B r]'-'H+l~ were handled in a similar way recently Q. A sizable gluonic content in the 77' 
meson has been indicated from a phenomenological analysis of the relevant data [6j . All these previous studies 
motivate us to make a more definite estimate of the gluonic contribution in the B — > ry^'^ form factors. 

In this paper we shall calculate the gluonic contribution to the B 77''^ form factors in the large-recoil region 
K/i \ using the perturbative QCD (PQCD) approach 0, H, [§]■ This approach is based on kr factorization theorem 
. [13, El J so the end-point singularities do not exist. It has been proposed to extract this gluonic contribution 
d ' from the measured B ij^'^lv decay spectra [12], which are, however, not available yet. To proceed with the 
calculation, we need to specify a scheme for the 77-77' mixing. After comparing the quark-flavor basis ^13j and the 
conventional singlet-octet basis, we adopt the former, in which fewer two-parton twist-3 rj^'"' meson distribution 
amplitudes are introduced. To reduce the theoretical uncertainty from the distribution amplitudes, we employ 
the Gegenbauer coefficients constrained by the data of both exclusive processes ^Sj, jM, 15], such as the r/^'^ 
transition form factors, and inclusive processes [IB], such as T(15') rj'X. It will be shown that the gluonic 
contribution is negligible in the B r] form factors, and reaches few percents in the B rj' ones. 

Whether the observed B rjK and B rj'K branching ratios [13] can be accommodated simultaneously still 
attracts a lot of attentions [H, [l^. We shall elaborate the impact of the gluonic contribution in the B rj^'^ 
form factors on this issue in the QCDF and PQCD frameworks. As noticed in [3], the gluonic contribution 
increases the branching ratios B{B il'K), but decreases B{B tiK). Since the QCDF predictions for 
both B{B — > r}' K) and B{B — > 77^) fall short compared with the data |^4], the gluonic contribution does not 
help. On the contrary, there is much room for this contribution to play in PQCD: the flavor-singlet amplitudes 
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were not taken into account in the earlier PQCD analysis of the B — > rj^'^K decays ^2011, whose predictions for 
B{B rj' K) [B{B tj^^)] s-re lower (higher) than the measured values. If stretching the gluonic contribution, 
it is likely to accommodate the data of the B rj^'^K branching ratios in PQCD. 

In Sec. II we compare the quark-flavor and singlet-octet schemes for the 77-77' mixing, and obtain the chiral 
enhancement scales associated with the two-parton twist-3 distribution amplitudes in both cases. In Sec. Ill 
we derive the factorization formulas for the quark and gluonic contributions in the B 77*-'^ form factors, and 
perform the numerical evaluation, together with a detailed uncertainty analysis. Our results are then compared 
with those obtained in the literature. The impact of the gluonic contribution on the accommodation of the 
measured B rj^'^^K branching ratios is discussed. Section IV is the conclusion. 



II. 77-77' MIXING AND DISTRIBUTION AMPLITUDES 

For the 77-77' mixing, the conventional singlet-octet basis and the quark-flavor basis have been proposed. 
In the latter the qq = {uu + dd)/^/2 and ss flavor states, labelled by the 77^ and 77^ mesons, respectively, are 
defined. The physical states rj and 77' are related to the flavor states through a single angle </), 

^ I'?) ^ - Trf^\ ( 1^?) 

with the matrix, 



fw= ^ (2) 

It has been postulated (Tpl that only two decay constants fq and need to be introduced: 

{0\srj5s\Vs{P)) = -^fsP^ , (3) 

for the light quark q — u or d. This postulation is based on the assumption that the intrinsic qq (ss) component 
is absent in the 77^ (77^) meson, ie., based on the OZI suppression rule. The decay constants associated with the 
77 and 77' mesons: 

(0|97''759|r/W(P)) = --!=/^%P^ 

{Q\-sri,sW'\p)) = (4) 

are then related to fq and fs via the same mixing matrix, 

( /'^ ^ ^ ^^''^^ f ^ 

Employing the equation of motion, 



f^, ) - ^^^^ I /. ) • 



dMl^'l^l) = 27777, qroq + ^ G^. G^'' , (6) 

and the one corresponding to the s quark, where G is the field-strength tensor and G the dual field-strength 
tensor, one derives the relation p^ . 

Ml - U\m^U{ct^) . (7) 

In the above expression the mass matrices are given by, 

ml 
ml 



m2 = f mlq + {V2/fq){0\asGG/{A7T)M (l//,)(0|a,GG/(4^)|77,) , 

«^ ^ (x/2//,)(0|a.GG/(4^)|77.) 777^ + (l//.)(0KGG/(47r)|77.^ ' ' ^' 
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with the abbreviations, 



rrigg = —{0\mumj5U + mddtj5d\r]q), 

Jq 

2 

mis = -r{0\mssi'^5s\vs) ■ (9) 

Is 

The above matrix elements define the chiral enhancement scales associated with the two-parton twist-3 rjq and 
r]s meson distribution amplitudes. Note that the axial U(l) anomaly |21,J is the only source of the non-diagonal 
elements of Mqs in the quark- flavor basis, also a consequence of the postulation that leads to Eq. ([3]). 

We repeat the above formalism for the 77-77' mixing in the singlet-octet basis, where the {uu + dd + ss)/\/3 
and {uu + dd — 2ss)/VQ states, labelled by the 771 and 773 mesons, respectively, are considered. From Eq. (P) 
and the quark content of the rji and rjg, mesons, we have the decompositions of the rj and 7/ meson states in the 
singlet-octet basis, 

''^ ^ = m ( \:i ) , (10) 



\v) J V 1^1 . 

with the angle = (p — 9i, 9i being the ideal mixing angle associated with the matrix, 

The decay constants defined in terms of the SU(3) flavor-singlet and flavor-octet axial-vector currents, 

(o|J^5h^'H^)> = -*/;(o^^ (12) 

for i = 1 or 8, are related to fq and /s through 
Compared to the conventional parametrization. 



with the matrix. 



Usi = ( ) ' (15) 

' smt^g costal ' ^ ' 



the decay constants of the 771 and 778 mesons, /i and /§, respectively, then connect to fq and fs 
Following the similar procedure, we derive the version of Eq. ([7]) in the singlet-octet basis, 



with the mass matrix. 



and the abbreviations. 



M|i = uH9)M^Usi , (16) 



+ {V3/fi){0\asGG/{47T)\m) , 
«i ^ ml, 777?i-f (V3//i)(0|a.GG/(4^)|77i) ' ' ^ 



77i|g = — — {0\muui^5U + ruddi^^d - 2msSi^5s\r]s) 
V6/8 
2 



2 



\muuij5U + viddi^bd -\- mssi^5s\7]s) , 
\muuij5u + mddij^d - 277^5 5775 s 1 771) 
771^^ = {0\muUij5U + ruddij^d + mssi^5s\7]i) ■ (18) 
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The above matrix elements define the chiral enhancement scales associated with the two-parton twist-3 iji and 
778 meson distribution amplitudes. 

Note that the four hadronic matrix elements m^q, rn^si {0\Q-sGG/ (47r)|77q), and (0|Q;sGG/(47r)|77s) are fixed by 
Eq. ([7]) in the quark-flavor basis, given the three inputs fq, fg, and (f>. However, contains six matrix elements, 
which can not be fixed completely by Eq. p6p . In fact, two OZI suppressed matrix elements {0\muUi"f5U + 
iTT-ddil5d\r]s) and {0\msSi^5s\riq) have been dropped in Il3l| . In the singlet-octet basis an approximation can be 
made to reduce the number of matrix elements [2^ l23j: the terms proportional to the light quark masses 77i„ 
and TTid are negligible compared with those proportional to rrig in Eq. (fT8|) . which then becomes 



mo 



V6/8 

2 



[0\msSij5S\r]s) I "^18 
\mssij5s\r]i) , m| 



^81 



V2/i 



/s 2 



(19) 



The three inputs fq, fs, and (j) in Eq. (O have been extracted from the data of the relevant exclusive processes 



-(1.07 ±0.02)/, , = (1.34±0.06)A, = 39.3° ± 1.0° 



which correspond, via Eqs. p3p and ([H]), to 0| 

/s « 1.26/, , /i « 1.17/, 



-21.2° 



-9.2° , 



(20) 



(21) 



in the singlet-octet basis. The above values are well consistent with those in 23], which were also extracted 
from the relevant data but based on the approximation in Eq. (jl9p . We stress that the approximations in the 
quark-flavor basis ^13i] and in the singlet-octet basis |23j for decreasing the number of hadronic matrix elements 
are very different. Therefore, the above agreement implies that these approximations make sense, and that the 
two bases will be equivalent to each other, if one chooses the parameters obeying the constraints in Eqs. ([7]), 
((131), (HID, (dH), and dUD for a calculation. 

Viewing Eqs. ^ and (fT8)l . it is obvious that fewer two-parton twist-3 distribution amplitudes are introduced 
in the quark-flavor scheme [5^. Hence, we adopt this scheme for the 77-77' mixing, and specify the 77^ and 
77s meson distribution amplitudes. Their two-parton quark components are defined via the nonlocal matrix 
elements, 



{r,qiPW^{z)q1,mO) 



(r;s(P)|.-(z)4(0)|0) 



eafc 



dxe"^-" {[75 /P]/3^0^(a;) + [75]/37™o<(a^) 



+^"o[75(A+ A- - l)]/37<^^(2;)} ' 



srab 



dxe'^P-'- {[75 /P]/37C(a;) + [l^h^m'^^'s {x) 



+^"o[75(A+ /i- - l)]/370f (a;)} , 



(22) 



where P = (P+,0, 0^) is the 77^^^ meson momentum, the light-like vector z ~ (0, z^,Ot) the coordinate of 
the q and s quarks, the dimensionless vector 714. — (1,0, Ot) parallel to P, 7i_ — (0,1, Ot) parallel to z, the 
superscripts a and b the color indices, the subscripts 7 and (3 the Dirac indices, Nc — 3 the number of colors, 
and X the momentum fraction carried by the q and s quarks. The chiral enhancement scales mg and ttiq have 
been fixed by Eq. ([7]), whose explicit expressions are 



mlq 


1 


2mq 


2mq 


™L 


1 






2ms 


2ms 



777,^ COS^ I 



2 -21 ^/^Js I 2 
771^, sm 0— '™ 



m.^) C0S(psm( 



777^, cos^ (j) + sin^ 1 



' (?Tt^^, — 777,^^ ) cos (j) sin (f> 



\/2/, 



(23) 



respectively, assuming the exact isospin symmetry 777.^ 
breaking effect, refer to [25j). 



md (For the inclusion of the isospin symmetry 
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The distribution amplitudes </)^^ are twist-2, and (f)^ ^ and (f>^ ^ twist-3. As explained in [26j, both the twist-2 
and twist-3 two-parton distribution amplitudes contribute at leading power in the analysis of exclusive B meson 
decays. We follow the parametrization for the pion distribution amplitudes proposed in [27l |. 



(/'tis) (x) 



2V2iVc 

2^/27v: 



1 + ( 30% - ^P^(,) ) Cy^2x - 1) 
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'?3^3 + ;3^P^(.)(l + 6af))|C;'^(2a;-l) 

1 7 n 3 , 



2J2N, 



:(1 - 2x) 



20' 
1 + 6 (5773 



20 



-Pq{s 



(1 - 10a: + lOx^) 



with the mass ratios pq = 2mq/mqq and ps = 2ms/mss, and the Gcgenbauer polynomials, 



The values of the constant parameters Oj, rj^ and uj^ in Eq. 



Cl^^it) = ^ (3-3Gi2 + 35i^ 



(24) 



(25) 



2, //a ciiiu 111 i^v^. will be given in the next section. 
The leading-twist gluonic distribution amplitudes of the rjq and rjs mesons are defined by [lB| 



{r,q{P)\Al{z)Al^iOm 



V2fq Cf 



sab 



V3 4V3A^2-1 



n_ • P 



dxe 



ixP-z 



fs Cf 



sab 



V34V3A^2-1 



nPP'' 



dxe 



ixP-z 



0^(^) 

x{l — x) 
x{l — x) ' 



with the notation A'^^{z)Al^{'w) = [Al{z)Al{w) - A'^{z)A''^{w)]/2 and the function 0,1281,12 

^^(.)(^) = ^'(1 - x)2i3|(^)Ci^/^(2x - 1) , Cl/\t) = 5t . 



mm 



(26) 



(27) 



According to the above definition, the gluon labelled by the subscript p. carries the fractional momentum xP. 
The two Gegenbauer coefficients i?| and -Bf could be different in principle. However, due to the large uncertainty 
in their values, it is acceptable to assume i3| = _B| = i?2. As shown later, the contribution from the above 
gluonic distribution amplitudes is smaller than that from the quark distribution amplitudes. Therefore, the 
subleading-twist gluonic distribution amplitudes of the 77^.3 meson will be dropped below. 

The gluonic distribution amplitude of the 77' meson defined in 133] is related to those in Eq (pS)) through 
Eq. HD, 



{v\Al{z)Al^mO) \ 
W\Al{z)Al^mO) I 



(r7,|^f^(z)A^](0)|0)\ 
{^s\Al{z)Al^{0)\0) 



(28) 



which also defines the gluonic distribution amplitude of the i] meson. Besides, our parametrization of the 
gluonic distribution amplitudes are identical to those in Q except a different definition of B2: their Gegenbauer 
coefficient is 2/9 times of ours. 



III. B rj'-'^ FORM FACTORS 

After defining the quark and gluonic distribution amplitudes of the rjq and rjs mesons, we are ready to calculate 
the B rj^'") transition form factors at leading order of the strong coupling constant as- In the B meson rest 
frame, we choose the B meson momentum Pi and the ry*-'^ meson momentum P2 in the light-cone coordinates: 

Pi = ^(1,1,0t), P2 = ^(p,0,Ot), (29) 
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where the energy fraction p carried by the ry^'^ meson is related to the lepton-pair momentum q = Pi — P2 
via = (1 — p)m?Q, thb being the B meson mass. The r/^'^ meson mass, appearing only in power-suppressed 
terms, has been neglected. The spectator momenta fci on the B meson side and k2 on the r?^'' meson side are 
parameterized as 



TUB 



0,a;i^=r,kiT , fc2= a;2p^=r,0,k; 



ruB 



t2T , 



Xi and X2 being the parton momentum fractions, and kiT and \<L2T the parton transverse momenta. 
We first compute the B Vq{s) form factors defined by the local matrix elements. 



(P1+P2) 



m; 



9m 



m. 



9m 



K(.)(P2)|6(0)ia'^''(?.d(0)|B(Pi)) 



TUB + 



"nQ(s) 



{m 



m. 



(30) 



(31) 



where the 77^(3) meson mass mn^^^^ will be set to zero eventually. The form factors F+_o are associated with the 

scmilcptonic decay B i]^'Hv, and Ft with B ri^'^l~^l~. For the involved b — > m, J transitions, the above 
form factors are decomposed into 



+(0,T) 



pBrig 



+ F, 



g+(0,T) ^ ^ 9+(0,T) 



(32) 



That is, the rjs meson state contributes only through the flavor-singlet pieces -F^gOgT- '^^^ ^ ~^ ^^'^ form 
factors are then obtained from the mixing, 



U{<t>) 



T^B-l)q 

+ (0,T) 
TpBris 
+ (0,T) 



(33) 



It is then expected that the gluonic contribution is more significant in the B rj' form factors than in the 
B ^ T] ones, since those from the rjq and r]s mesons add up in the former, but partially cancel in the latter. 



B 



(') 



FIG. 1: Gluonic contribution to the B ■q'^'^ form factors. Another diagram with the two gluons crossed is suppressed. 



The factorization formulas for ^"^''^0 qT similar to those for the B n form factors [1^ [32 



F^^''{q^) = — 7rm|C_F / dxidx2 / bidbib2db2<^Bixi,bi) 



(1 + X2p)(t)q{x2) +rq[--l-2x2]rq ix2) + rq{l - 2x2)(t)g {X2) 



E{t^^'^)h{xi,X2MM) 



+2rg^PE{t^^)h{x2,Xi,b2,b,) \ , 



Pqo''{<l^) -^^itT'bC^fP J dxidx2 J bidbib2db2(j)B{xi,bi) 



(1 + X2p)4'g{x2) + Tq (1 - 2x2) (t)g {X2) + [ - - I ^ 2x2 ] 4>g ix2) 



E{t^^'^)h{xi,X2MM) 



+2rgC^PE{t^^)h{x2,Xi,b2M)^ , 
F^T^{q'^) = -^niTLgCF J dxidx2 J bidbib2db2(j)B{xi,bi) 



(X2) [--+X2] (t)q {X2) - rgX2(f>g {X2) 



E{t'^^^)h(xi,X2MM) 



+2rg^^E{t^^)h{x2,xi,b2,bi) } , 



(34) 



with the color factor Cp — 4/3, the B meson wave function (j)B, the impact parameter bi (62) conjugate to the 
parton transverse momentum kix (A;2t), the mass ratio rg ~ rriQ/mB, the hard function h, and the evolution 
factor, 



(35) 



The choice of the hard scale t, and the explicit expressions for h, for the Sudakov exponent Sb associated 
with the B meson, and for the Sudakov exponent Sr/^ associated with a light meson are referred to [26]. The 
threshold resummation factor associated with the hard function is the same as in the i? — > tt form factors j34| . 
The coefficient l/\/2 appears, because only the u or d quark component of the rjg meson is involved. We point 
out that the term proportional to 2/ p in F^r^'' is missed in [3^ . 

For the gluonic contribution, it has been argued Q that the diagram in Fig. [T] with the two gluon emitted 
from the light spectator quark of the B meson is leading. Another diagram with the two gluons crossed, giving 
the identical contribution, is not displayed. The third diagram vanishes, in which the virtual gluon from the 
u and u quark annihilation couples the two valence gluons in the jy*^'-* meson. The flavor-singlet pieces in the 
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B ^ rjq form factors are written as 



xxi[l + (1 - p)x2)]E{t(^)h{x2,xi,b2,bi) , 
FgQ '{q ) = -g7rmg/g _ i P J dxidx2 j 6id6i62d620B(a;i, 61)— ^ 



.(1-2^2) 

xxi[l - (1 - p)a;2)]£'(t^^)/i(a;2,a;i,62,&i) , 

^gT il ) = -gTrmg/^ jV^ - 1 y "^i"^2 y bidbib2db2(pB[Xi,bi) ^ {I - x ) 

xxi{l + X2)E{t^^)h{x2,xi,b2,bi) . (36) 

A factor 2 has been included, which is attributed to the identical contribution from the second diagram men- 
tioned above. The calculation is similar to that of the r/'g*g^*^ vertex in |36, 37]. The expressions for ^^''^q gT 

are the same as in Eq. but with fq being replaced by fs/V2. It is observed that Eq. is proportional to 
the small momentum fraction xi ^ A/ms '38|, A being a hadronic scale, compared to the quark contribution. 
Because the above factorization formulas do not develop the end-point singularities, if removing fcy, the gluonic 
contribution can in fact be computed in coUinear factorization theorem. 
The evolution factor in Eq. (j36p is given by 

E{t) = asit)e-^^^'^-^^^'^ , (37) 

where the Sudakov exponent 5*^ is associated with the gluonic distribution amplitudes of the ?7g.s mesons. 
Following the studies in [s^, its expression is, up to the leading- logarithm accuracy, similar to 5^,^, but with 
the anomalous dimension agCp/'!^ being replace by asC^i/Tr, Ca — 3 being a color factor: 



Sait) = SG{x2P^,b2) + SG{a~X2)P^,b2) 



1/b ^^ 



In ( ^ ) A{as{^Ji)) 



A - ^Ca ■ (38) 



That is, the Sudakov suppression is stronger in the gluonic distribution amplitudes than in the quark ones. There 
is no point to include the next-to-leading-logarithm resummation, since the gluonic distribution amplitudes 
are not very certain yet. For consistency, we neglect the single- logarithm renormalization-group summation 
governed by the anomalous dimension of the gluon wave function. We adopt the one-loop expression of the 
running coupling constant as, when evaluating the above Sudakov factors. 

We then perform a detailed numerical analysis, including that of theoretical uncertainties fioj. The B meson 
wave function is the same as in .4J^ with the shape parameter varied between cub = (0.40 ± 0.04) GeV. There is 
another B meson wave function in the heavy-quark limit , whose contribution to transition form factors may 
be finite. However, it has been shown that its effect can be well mimicked by a single B meson wave function, 
if a suitable ojb is chosen [i^l- Therefore, we adopt this approximation, and vary lob in the above range. To 
obtain the chiral enhancement scale TOq, we need the masses to^ — 0.548 GeV and to^/ = 0.958 GeV, and 
the inputs in Eq. ([^0]) . Because meson distribution amplitudes are defined at 1 GeV, we take the light quark 
mass mq{l GeV) = (5.6 ± 1.6) MeV [i^. We have confirmed that the twist-2 and twist-3 quark distribution 
amplitudes of the 77' meson in [33| are consistent with Eqs. p6|) and (|22p . if their 77' meson decay constant /,,' 
is regarded as f^, in the singlet-octet basis. Hence, it is legitimate to adopt = —0.008 ±0.054 extracted from 
the relevant data (l6j for the twist-2 distribution amplitude. The twist-3 distribution amplitudes have not yet 
been constrained experimentally, so we choose 773 = 0.015 and lu^ = —3 the same as for the pion distribution 
amplitudes [2^ . It is not necessary to specify the parameters involved in the quark distribution amplitudes of 
the 77s meson here. The overall coefficients in Eq. (1^51) have been arranged in the way that the decay constants 
in [121 and in Eq. satisfy Eq. (US]). Following Eq. ((281), the range of B2 = 4.6 ± 2.5 extracted in [H] can 
also be adopted directly. The theoretical uncertainty arising from the variation of the above parameters will be 
investigated. 
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^ + ,0 


pBr, 






m 


0.147 


0.139 


0.121 


0.114 


ratio 


0.0031 


0.0028 


0.023 


0.021 


Fi 


0.252 


0.237 


0.252 


0.237 


F2 


0.0034 


0.0029 


0.0034 


0.0029 



TABLE I: Form factor values at maximal recoil, the ratios of their gluonic contributions, and values of _Fi,2 defined in 
Eq dSl}. 



The following parametrization for the B rj'^'^ form factors was proposed in 



(39) 



where Fi {F2) corresponds to the quark (gluonic) contribution, and a factor 1/V2 is included to match our 
convention. Fi has been set to the B —* tt form factor, and the unknown F2 varied arbitrarily between [0, 0.1] 
The above parametrization also applies to the other form factors and Ft- It is easy to identify, from 
Eq. ([33]) . the relations of Fi,2 to our form factors. 



V2U 



Ff''(0), F2 



V3U 



(40) 



To have a picture of the magnitude of the gluonic contribution, we present in Table U the central values of the 
form factors at maximal recoil, the ratios of their gluonic contributions, and the values of Fi 2- It is indicated 
that the gluonic contribution is about 0.3% in the B rj form factors, and about 2% in the B rf ones. The 

central value of Fi for the form factors F^^ is indeed close to the B tt form factors F^q{0) [2^. The central 

value of F2 = 0.0034 for ' is located within [0, 0.1] [4], but near the lower end of the interval (to be precise, 
F2 runs between [0.0016,0.0052], considering the range of the parameter ^2). We notice the difference between 

F2 for F^^ (0) and F2 for F,^^ (0), which is obvious from the factorization formulas in Eq. ([5^ : -F^'''(0) 

contains an additional term proportional to X2 compared with f^gQ(O). Our observation is thus contrary to 

the assumption F^''_^q(0) — F^^'" {Q) postulated in the analysis of the B —t r]'^'h^l^ decays [S]. Nevertheless, 
this difference is not crucial for the estimate in Q, because the gluonic contribution is not dominant. 




(a) 



10 

(GeV^) 



0.5 
0.4 
0.3 
0.2 
0.1 




(b) 



10 

(GeV^) 



J, and (b) -F^g ti corresponding to the central values of the inputs with the solid. 



FIG. 2: (j^ dependence of (a) F 

dotted and dash-dotted lines representing the F+, Fq, and Ft , respectively. 



The dependence of the B ^ rj and B —f ?/ form factors corresponding to the central values of the inputs 
are displayed in Figs. [2 (a) and (b). The ratios of the gluonic contributions to the total values of the form factors 
are shown in Fig. O It is found that the gluonic contribution remains negligible in the B rj form factors, and 
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FIG. 3: dependence of the ratio of the gluonic contribution to the total values of F^^^^ (solid lines), F^^^-^ ' (dotted 
lines), and f!^^'^^ ^ (dash-dotted lines) corresponding to the central values of the inputs. 



about 2% in the B ^ rj' ones in the whole large- recoil region. The form factors -Fjf''' \q'^) exhibit a smaller 
slope with q^, because the overall coefficients of F^^^q contain the energy fraction p as shown in Eqs. ([M]) and 
consistent with the large-energy form- factor relation in The B ^ rj form factors are proportional to 
cos 4>Fq ^ and the B —>■ rj' form factors to svcKpFq^'' plus some amount of gluonic contributions. Because of 
the small gluonic contribution indicated in Fig. [31 we have F^^{q^) > F^^ {q^) shown in Fig. [2] simply due to 
cos^ > saKp for (j) « 39.3°. This observation is in agreement with the tendency exhibited in the measurement 
of the semileptonic branching ratios [ist . 

B{B+ r^l+n) = (0.84 ± 0.27 ± 0.21) X lO"'^ < 1.4 X 10""* (90% C.L.) , 
B{B+ -> i^'l+vi) = (0.33 ± 0.60 ± 0.30) X 10""^ < 1.3 X 10""* (90% C.L.) . (41) 

More precise data will provide the information on the importance of the gluonic contribution. Our analysis 
implies that the of dep endence of the gluonic contribution is weaker than that of the quark contribution. Hence, 
the assumption [5|, Il2l | that both pieces show the same dependence is not appropriate. 




FIG. 4: Dependence of ttIq (a) on (f> in units of degrees, and (b) on fs/ fq- 



We then investigate the theoretical uncertainty for the B ry*^'-' form factors from each of the input pa- 
rameters, and the results are presented in Table HTl The form factors F^^''^^ ^j, in Eq. (p4)) receive substantial 
contributions from the twist-3 distribution amplitudes (j)^ and (j)^ , which appear together with the mass ratio 
rq = 171^1 niB- The chiral enhancement scale TOq changes rapidly with the mixing angle <f> and with the decay 
constants fq^s as illustrated in Fig.[4l This explains the sensitivity of the form factors to these inputs. We notice 
that mg runs into the unrealistic negative region easily, as increasing fs, leading to the negative form factor 
values in Table |TT1 Therefore, we prefer not to vary fs- The dependence of the form factors on the Gegenbauer 
coefficients a\ and B2 is weak, because a\ is close to zero and the gluonic contribution is subdominant. The 
variation of the ratios of the gluonic contributions in the B r]^'^ form factors with each of the parameters 
is listed in Table IIIII Ignoring the effect from changing fs for the reason stated above, we conclude that the 
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UJB 


= (0.36 - 0.44) 


0.174 


- 0.127 


0.164 


- 0.119 


0.127 - 0.106 


0.136 


- 0.099 


h-- 


= (1.05 - 1.09)/. 


0.057 


- 0.239 


0.049 


- 0.230 


0.049 - 0.198 


0.042 


- 0.190 


fs -- 


= (1.28 - 1.40)/. 


0.359 - 


(-0.063) 0.348 - 


(-0.070) 0.296 - (-0.049) 0.287 - 


(-0.055) 


<^ = 


= (38.3 - 40.3) 


0.093 


- 0.208 


0.084 


- 0.199 


0.076 - 0.179 


0.069 


- 0.171 


niq 


= (4.0 - 7.2) 


0.191 


- 0.124 


0.181 


- 0.116 


0.158 - 0.104 


0.150 


- 0.097 


al 


= (-0.062) - 0.046 


0.144 


- 0.152 


0.136 


- 0.143 


0.121 - 0.126 


0.113 


- 0.119 


B2 


= 2.1-7.1 


0.148 


- 0.148 


0.139 


- 0.139 


0.122 - 0.125 


0.115 


- 0.117 



TABLE II: Theoretical uncertainty of the form factors at maximal recoil from each of the parameters. 



gluonic contribution is negligible, always below 1%, in the B ^ rj transitions, and may reach order of 10% in 
B rj'. 
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0.29 - 0.34 
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2.2 




= (1.05 - 1.09)/. 


0.75 - 0.21 


0.73 - 0.18 


5.7 - 


1.5 


5.6 - 


1.3 


fs-- 


= (1.28 - 1.40)/. 


0.14 - (-0.64) 0.12 - (-0.49) 0.94 - 


(-6.0) 0.81 - 


(-4.5) 


= 


(38.3 - 40.3) 


0.55 - 0.20 


0.51 - 0.17 


3.7 - 


1.6 


3.5 - 


1.4 




= (4.0 - 7.2) 


0.24 - 0.37 


0.21 - 0.33 


1.8 - 


2.7 


1.6 - 


2.5 


al 


= (-0.062) - 0.046 


0.32 - 0.30 


0.29 - 0.27 


2.4 - 


2.2 


2.1 - 


2.0 


B2 


= 2.1 - 7.1 


0.14 - 0.48 


0.13 - 0.43 


1.1 - 


3.5 


0.95 


- 3.1 



TABLE III: Theoretical uncertainty of the ratios (in %) of the gluonic contributions in the form factors at maximal recoil 
from each of the parameters. 



At last, we discuss the impact of the gluonic contribution in the B ij'^'^ form factors on the predictions 
for the B rj^'^K branching ratios in QCDF and in PQCD. The current data of the branching ratios are 
summarized below [l7j : 

B{B^ ^ ry'ifi) = (69.7l|j) x 10-*^ , 
B{B°->7fK°) = (64.9 ± 3.5) X 10~^ , 
B{B^^r]K^) = (2.2±0.3) X 10~\ 

rjK") < 1.9 X 10"^ . (42) 

Table 4 in the QCDF analysis [3] shows the dependence of the B —> rj^'^K branching ratios on the gluonic 
contribution F2: B{B — > v'^) increases with F2, but B{B rjK) decreases. The reason is as follows. The 
gluonic contribution enhances the B —f rjqK amplitude (containing the B rjq transition), such that its 
cancellation with the B — > Krjs amplitude (containing the B K transition) becomes more exact in the 
B rjK decays. However, the above two amplitudes are constructive in the B r/K decays, whose branching 
ratios then exhibit an opposite behavior with the gluonic contribution. Table 4 in [J| also shows the predictions 
from the default scenario of inputs (with the strange quark mass nis — 100 MeV), B{B^ v'^^) ~ 42 x 10~^ 
and B{B^ V^^) ~ 1-7 x 10^® for F2 — 0, both of which fall short compared to the data. Enlarging F2, the 
predicted B{B^ -q' K^) increases, but B{B^ V^^) decreases and deviates more from the measured value. 
That is, there is no much room for the gluonic contribution to play in QCDF. Nevertheless, a smaller = 80 
MeV does help, since it Hfts both B{B'^ — > -q' K^) and B{B^ -qK^) as demonstrated in Table 4 of |4i]. 

We have stated that the flavor-singlet amplitudes were not taken into account in the PQCD study of the 
B rj^'^K decays The predictions B{B° r]'K°) 45 x lO'^ and B{B° r]K°) « 4.6 x 10"^ were 

obtained for rris = 100 MeV and for the chiral enhancement scale ttiq = 1.4 GeV (see Table 2 of [l^l), which is 
within our parameter range for TOq as displayed in Fig. |4l Because of the dynamical enhancement of penguin 
contributions [1, E^, both the above branching ratios are larger than those in QCDF from the default scenario 
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[J]. Comparing the PQCD predictions with the data in Eq. (|^^ . and knowing their dependence on F2, there is 
more room for the giuonic contribution to play apparently. Though the central value of the gluonic contribution 
is small, we can not exclude the possibility of accommodating the observed B r]^'^K branching ratios in PQCD 
under the current theoretical uncertainty. The complete PQCD analysis of the B rj^'^K decays including all 
flavor-singlet amplitudes will be published elsewhere. 

Before concluding, we mention that the form factor F^^(0) = 0.16 ± 0.03 has been derived from light- 
cone QCD sum rules 47[, which, however, did not include the two-parton twist-3 and gluonic contributions. 
Instead, the three-parton rj meson distribution amplitudes were taken into account. Even so, their result is 
in agreement with ours in Tables U and |TT1 The form factors F^"^" (0) = -0.023 ± 0.048 (0.045 ± 0.086) and 
i^^''»(0) = -0.099 ± 0.024 (-0.066 ± 0.043) [48| have been extracted by fitting parametrization in the soft- 
coUinear effective theory p9l | to the data of two-body nonleptonic B meson decays, where the numbers in the 
parentheses represent the second solution of the fitting. These values, after considering the large uncertainties, 
still differ from our observations dramatically, F^'''" (0) = 0.190 and Ff''=(0) = 0.0005. Our Ff''=(0), containing 
only the gluonic contribution at leading order of as, is much smaller than _F^''''(0). 



IV. CONCLUSION 



In this paper we have calculated the gluonic contribution to the B — > 77*^'^ transition form factors in the large- 
recoil region using the PQCD approach. The quark-fiavor and singiet-octet schemes for the 77-77' mixing were 
compared, and it was found that fewer two-parton twist-3 distribution amplitudes could be introduced in the 
former. The leading-twist quark and gluonic distribution amplitudes of the rig and 77^ mesons in the quark-flavor 
basis were constrained experimentally. The parameters involved in the two-parton twist-3 quark distribution 
amplitudes were determined either by equations of motion associated with the 77-77' mixing, or taken the same 
as in the pion distribution amplitudes from QCD sum rules. Therefore, we are able to predict the unknown 
gluonic contribution with less theoretical uncertainty. It has been shown that this contribution is negligible in 
the B —f r] form factors, and reaches few percents in the B r]' ones. These predictions can be confronted 
with the future measurement of the B rj^'^lv decay spectra. We have elaborated the impact of the gluonic 
contribution on the B — > rj^'^K branching ratios obtained in QCDF and in PQCD. The observation is that the 
gluonic contribution does not help accommodating the measured B — s- rj^'^K branching ratios in QCDF, but 
does in PQCD. 
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